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PURPOSE: To examine and compute the distortion experienced by a wide-band frequency- 
modulated (or phase-modulated) signal as it is transmitted through the earth's 
ionosphere. The results can be used as an aid in designing satellite systems that 
transmit television or multichannel telephone signals. 

RELATED TO : RAND's continuing study of communication satellite technology for the 
National Aeronautics and Space Administration. 

DISCUSSION AND METHODOLOGY : When a terrestrial station is linked to an earth satel- 
lite at synchronous altitude, the propagation path traverses virtually the entire 
ionosphere in a direction depending on the relative location of the station and 
the satellite. This path can range from the zenith to the near horizon. The aver- 
age integrated electron density through the ionosphere is about 10^ electrons/m^, 
but can vary by a factor of 10. Efficient use of satellite power requires wide- 
band transmissions, and distortion of signal occurs because of the dispersive ef- 
fects of the ionosphere, with greater effects at lower carrier frequencies and 
during periods of solar activity. Previous work has been concerned with the trans- 
mission of individual pulses, and with spectral analyses for wide-band FM, based 
on expanding the quasi-stationary approximation to the filtering effect of the 
ionosphere and retaining the lower-order distortion terms. The present analysis 
used a modification of the Bedrosian-Rice technique (P-3502, Distortion and Cross- 
talk of Linearly Filtered, Angle-modulated Signals , March 1967). An expansion 
was obtained for the phase of an angle-modulated (AM) signal passed through a 
general linear filter, and a spectral analysis was performed. The ionosphere was 
then modeled as such a filter, and numerical computations of the output linear- 
signal, cross-power, and intermodulation spectra were performed for the case of 
an FM signal having a flat gaussian baseband. Two examples of current interest 
were examined numerically: a high-capacity communication link of the Intelsat 
type, and a wide-band video transmission near the upper end of the UHF television 
band. 

PRINCIPAL FINDINGS : Ionospheric dispersion will not cause significant intermodulation 
distortion in a typical wide-band FM communication satellite link except under the 
most severe conditions. Nevertheless, for a given baseband signal and r-f band- 
width, operation over a given path should be at as high a carrier frequency as 
possible. 

The situation is less promising with respect to color television. The calcu- 
lated ratios of signal-to-distortion and to cross-talk are about 20 db. It is dif- 
ficult to assess the subjective effect of intermodulation distortion for a hypo- 
thetical TV distribution system using wide-band FM transmission from a satellite 
relay as compared with the "snow" produced by a conventional vestigial-sideband AM 
transmission corrupted by thermal noise. The thermal noise in AM video is uniform, 
whereas the intermodulation distortion of FM varies roughly as the square of the 
baseband frequency and can display a strong correlation with the signal. Color 
transmission would suffer also from a degradation of color separation because of 
synchronizing phase errors. Just how objectionable these effects would be cannot 
be determined by computation. An experimental investigation would be desirable. 

MW 
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PREFACE 

This study, undertaken as part of RAND' s continuing study of 
Communications Satellite Technology for the National Aeronautics and 
Space Administration, considers the intermodulation, or distortion, 
experienced by a wide-band, frequency-modulated signal as it is trans- 
mitted through the earth 1 s ionosphere. 

The efficient use of the power available from satellites requires 
the use of wide-band transmission. Distortion occurs because of the 
dispersive effects of the earth's ionosphere, with the effects greater 
at the lower carrier frequencies and during periods of solar activity. 
This Memorandum studies this phenomenon and provides a means for com- 
puting the resulting distortion. The results of the study can be 
used as an aid for designing satellite systems that transmit television 
or multichannel telephone signals. 
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SUMMARY 

Communication links to earth satellites are not usually affected 
adversely in traversing the ionosphere. However, the nonlinear vari- 
ation of the refractive index with frequency causes a dispersion which 
can result in significant signal distortion under certain circumstances. 
A wide-band signal operating at the lower carrier frequencies and taking 
a low-elevation path through a disturbed ionosphere constitutes such 
a case. The effect of ionospheric dispersion is examined in this Memo- 
randum for the particular case of frequency (or phase) modulation. 

In a previous analysis (P-3502, Distortion and Crosstalk of Line - 
arly Filtered, Angle-modulated Signals , The RAND Corporation, by 
E. Bedrosian and S. 0. Rice) an" expansion was developed for the output 
phase of an angle-modulated signal which had been passed through a 
generalized linear filter. The ionosphere is modeled as such a filter 
in this Memorandum, and a spectral analysis of the demodulated signal 
is performed assuming a gaussian modulating waveform. From this, 
formulas are obtained for the output signal-to-distortion and signal-to- 
cross-talk ratios. All of the spectral results are presented graphi- 
cally for FM with a uniform baseband. 

Two cases of current interest are examined numerically to illus- 
trate the use of the results. The first case typifies a high-capacity 
communication link of the Intelsat variety, and it is seen that a 
significant degradation can occur only under the most unusual condi- 
tions. The second example considers a wide-band TV transmission near 
the upper end of the UHF TV band; in this case, ratios of signal-to- 
distortion and to cross-talk of about 20 db result. Ratios of this 
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order are clearly cause for concern, but it is not possible to assess 
their subjective effect by relating them to the familiar interference 
produced by thermal noise which is uncorrelated with the signal. 
Also, a TV signal, particularly if color, is not as well modeled by 
a gaussian signal as is a typical multichannel telephone baseband. 
For these reasons, it is concluded that an experimental investigation 
would be desirable both to verify the theory and to develop subjective 
interference criteria. 
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I, INTRODUCTION 

Communication links between terrestrial stations and earth satel- 
lites or space probes invariably employ operating frequencies well above 
those for which ionospheric absorption or refraction are significant. 
Other effects, such as Faraday rotation and path splitting, are either 
found to be unimportant or are accommodated without difficulty. However, 
this is not necessarily the case for the dispersion produced by the non- 
linear variation of the velocity of propagation with frequency. A wide- 
band signal traversing a disturbed ionosphere can experience considerable 
distortion at the lower carrier frequencies. 

Much of the work to date on ionospheric-induced distortion has been 
concerned with the transmission of individual pulses or of pulse-type 
modulation, although some attention has also been directed to analog- 
type modulation. This work is briefly reviewed in Ref. 1. The earliest 

work which appears to relate directly to the results presented here was 
(2) 

done by Rice^ y in 1957. He approximates the expression for the output 
phase of a filtered signal for the case of small distortion and then 
performs a spectral analysis using second-order modulation terms. In 
particular, he obtains the spectrum of the "interchannel interf erence" 
(Eq. (5.5)) when an FM signal having a flat gaussian baseband is passed 
through a filter with a uniform attenuation and quadratic phase charac- 
teristic. It will be shown later in this study that this yields a narrow- 
band approximation to the more precise result presented here. 

Direct applications to ionospheric propagation of wide-band FM, 
which is of interest in broadcasting and relaying TV signals, were per- 
formed by Shaft, ^ and Prosin^ in 1963, and by Denton, Scheibe and 



Huntsberger v in 1965 using the quasistationary approximation to the 
filtering effect of the ionosphere. In each case, the distortion spec- 
trum (Shaft, Eq. (21); Prosin, Eq. (39); Denton, et al. , Eq. (4.26)) is 
obtained by expanding the quasistationary approximation and retaining 
the lower order distortion terms, thereby yielding results which agree 
with Rice's approximation given in Ref. 2. 

The following analysis is based on the technique described by 
Bedrosian and Rice^ ' wherein an expansion is obtained for the phase of 
an angle-modulated signal passed through a general linear filter. The 
expansion is then specialized to the case of a symmetrical band-pass 
filter, and the output spectrum is computed numerically for the case 
when an EM signal, having a flat gaussian baseband, is passed through 
a single-pole filter. 

In this Memorandum, a spectral analysis is performed when a general 
filter is used. The ionosphere is then modeled as such a filter and 
numerical computations of the output spectrum are again performed for 
the case of an FM signal having a flat gaussian baseband. The results 
are interpreted for specific examples of transionospheric propagation of 
a wide-band FM signal using communication and video basebands. The 
signal-to-distortion ratios predicted in a TV case of potential interest 
are shown to be poor enough in comparison with the accepted signal-to- 
noise criteria for high-fidelity transmission to warrant an experimental 
investigation to determine the nature and subjective effect of distortion. 
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II. OUTPUT PHASE FUNCTION 

Consider the angle-modulated signal 

s(t) = exp i[2TTf Q t + cp(t)] (1) 

where f is the carrier frequency and cp is the modulating signal. 

When this signal is passed through a filter having an impulse response, 

g, and a transfer function, G, where 



G(f) = 



. . -i2TTft , . . 

g(t)e dt, g(t) 



G(f)e i2nft df (2) 



are a Fourier pair, then the output phase, 9, is given by^^ 

CO 

9(t) = -0 + Re $(t) + V -V Im i n f (3) 
o n • n 

n=2 

where (3^ is the phase shift at the carrier frequency 

i3 Q = -Im log G(f ) (4) 
and the function § is given by 



4(t) = 



Y(T)9(t-T)dT (5) 



-i2TTf t 

v(0 - P(t) o(f o ) ° 

The first four coefficients in the series are 

2 

f 2 " V f 4 = F 4 " 3 F 2 

f 3 = F 3 , f 5 = F 5 - 10 F 3 F 2 (7) 



where 
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F = 
n 



Y(T)[cp(t-T)-$(t)] n dT (8) 



0 

Additional coefficients and a discussion of the convergence of this 
series are presented in Ref. 6. 

The impulse response, y> given by Eq. (6) is associated with a 
filter which is a normalized, frequency- shifted version of the original 
filter. The transfer function of this hypothetical filter is 

G(f+f ) 
o 

and it follows, from Eqs. (2) and (6), that 

y(T)dT = T(0) = 1 (10) 
0 

If the original filter is physically realizable, then its impulse 
response, g, vanishes for t < 0 and is real, while its transfer function, 
G, has an even real part and an odd imaginary part by virtue of Eq. (2). 
From Eq. (6), it is seen that the impulse response, Y, of the normalized 
frequency- shifted filter also vanishes for t < 0, but y is complex in 
general, and its associated transfer function, V , lacks the symmetry 
properties of G. 

It is convenient, therefore, to write 

Y(t) = Y r (t) + i YjU) (11) 

where, from Eq. (10), 



Y r (t)dt = 1, 



Y i (t)dt = 0 (12) 
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Then, since y^ and y are real functions, it follows that they can be 
associated with the transfer functions 



A(f) = 



. . -i2rrft J 
Y r (t)e dt 



A(f) = J Y^^e" 12 "" dt (13) 
0 

which do have the even real-part and odd imaginary-part symmetry of G. 
Of course, A and A are complex, in general. In terms of the real 
and imaginary parts of F, 

F(f) = r r (f) +.iF.(f) = A(f) + iA(f) (14) 

from which it is seen that 

A(f) = j[r r (f)+r r (-f)] + i|[r.(f)-F.(-f)] 

A(f) = \ [r.(f) + F.(-f)] - i 1 [r r (f)-F r (-f)] (15) 

These functions completely characterize the normalized, frequency-shifted 
filter to be used in the ensuing spectral analysis. 

The output phase through terms of third order is, from Eq. (3), 

9(t) = -f> Q + Re $(t) + jr Im i 2 f £ + Im i 3 f 3 + ... (16) 
Writing the function $ in terms of its real and imaginary parts yields 

$(t) = $ (t) + i §.(t) 



= J Y r (T)<P(t>T)dT + i J Yi (T)Cp(t-T)dT (17) 

o o 1 
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Substituting in Eqs. (7) and (8) and noting Eqs. (12) then leads to 



9(t) = -3 o + $ r (t) - -| 



( T ) 



^(t-T)-^^ (t) J dT 



V r (T) 



2 $ i (t) 



0 



J 

=F(t-T)-$ r (t) J dT 
Y i (T) [«P(t-T)-§ r (t)] dT + 



(18) 



which is the desired form of the output phase. The leading term 
yields the basic phase shift at the carrier frequency. Depending on 
the nature of the input phase, <p, additional phase shifts will, in 
general, be contributed by the remaining terms. Even if 9 is zero 
mean, phase shifts can be generated by terms such as the third one. 

The second term in Eq. (18) is a linearly filtered version of 
the input and represents the principal time- vary ing component of the 
output when the distortion is small. Successive terms yield the 
various orders of distortion to the output. 
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III, SPECTRAL ANALYSIS 

The spectral density, Wq , of the output phase, 0, will be deter- 
mined by the conventional technique of taking the Fourier transform of 
the autocorrelation function of 9 

W e (f) = 3 E 9(t)9(t-T) (19) 

where 3 denotes the Fourier transform operator and E, the expectation 
operator, denotes an ensemble average defining the autocorrelation 
function 

Rq(t) = E e(t)9(t+T) (20) 
and where R and W are a Fourier pair. In terms of phase rates 

W£(f) = (2TTf) 2 w Q (f); vyf) = W^(f)/(2Hf) 2 (21) 

If angles are measured in radians and frequency in hertz, the units of 

2 2 
Wq are rad /Hz, and those of Wq are (rad/s) /Hz. 

When the expansion for the output phase given by Eq. (18) is sub- 
stituted in Eq. (19), the results contain a variety of functional pair- 
ings. Let £ , §, f^, f^, etc., denote the terms in the output phase 
given by Eq. (18) and denote the individual pairings by terms such as 

§Xf = 3 E $(t)f (t+T) (22) 

Then, as discussed in Ref. 6, the 11 linear- signal 11 component of the out- 
put spectrum is identified as the §X$ term and denoted by W^. The $xf 
and fxf terms contribute some output spectral components in which the 
input spectrum, W , appears as a multiplier; these are identified as 
the " cross-power" component of the output spec t rum and are denoted by 

Q 

Wq • The " intermodulation 11 component of the output spectrum is denoted 
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by Wq and is given by the balance of the frequency-dependent spectral 
terms which result from the fXf pairings. In addition to these three 
types, which are identified in Ref. 6, a zero- frequency or "dc" term, 
denoted by 9 , appears in the analysis of the general case. 

For the purposes of the spectral analysis, the input signal, <p, 
is taken as a gaussian process having a spectral density, (the 
details of the analysis are presented in the Appendix). The leading 
terms of the various components of the output spectral density are 
given by 



2 2 

r = 3 - 23 

o o 



J dp W (p)Re A(p)A(-p) 



w£(f) = |A(f)| Z W 9 (f) 



w e (f) = % (f) 



dp W^(p){Re A(f)[A(p)A(-p-f)-A(p)A(-p-f)] 



+ |A(f)| 2 [|A(p)| 2 -|A(p)|] 2 + 2Re A(f)A(-f)A(p)A(-p)] 



W 9 (f) = \ I dP w C p(P) w cp ( f -P) l A (P) A ( f -P) + A(p)A(f-p) - A(f)| 2 (23) 



The principal dc term results from the 3 xj3 operation with a second- 
order correction from the 3*f£ term; third-order corrections, which 
are not shown, come from the 3 x f^ and ^2*^2 terms * As statec * above, 
the linear-signal component is due to the $X$ term. The cross-power 
component is due solely to the * x f^ term since the f 2 Xf 2 tem y ields 
only an intermodulation component in this case. 
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IV. APPLICATION TO IONOSPHERE 

A transfer function which is characteristic of the transionospheric 
propagation path is required to use the foregoing analysis. Since, for 
all practical purposes, the ionosphere behaves as a linear medium and is 
time- invariant on time scales of interest for communication purposes, it 
follows that the propagation path under consideration can be modeled as a 
cascade of the same path in vacuo and a filter which reflects only the 
properties of the ionosphere. The transfer function of that filter is 
then given by the ratio of the steady- state frequency response of the 
actual transionospheric path to that of the free-space path. 

The relative bandwidth of the transmitted signal will be assumed to 
be sufficiently small , and the carrier frequency sufficiently large in 
comparison with the plasma frequency, so that the effect of ionospheric 
absorption can be neglected. Since the distance dependence is the same 
for both paths, it follows that the transfer function depends solely on 
the difference between the phase characteristics of the two paths, and 
is given by 

G(f) = expCik J (n-l)ds] (24) 

where 

k = Znf/c is the phase constant 
g 

c = 3 x 10 m/s is the speed of light 

f 2 

n = ll-f /f j is the refractive index of the ionosphere at frequen- 
^ cies well above the plasma frequency, f , and the 

integration is along the propagation path. 

For frequencies large in comparison with the plasma frequency , the 

refractive index can be expanded in a binomial series. Then, expanding 

the leading term, l/f, in a Taylor's series about the carrier frequency, 

f , yields the high-frequency narrow-band approximation 
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0(f) - e*p{i 1 [f - + 1^] J f * ds } (25) 

of f r 

o o 

where the various terms represent, respectively, the constant, linear 
and quadratic components of the frequency-dependent phase shift intro- 
duced by the ionosphere. To simplify the expression, let 

2 C ^ 

c nff 2 ds- 
P 

Then, denoting normalization to f by underlining yields 

^o " " Im l0 § G < f o> " ^ ' fo = f (27) 

for the carrier-frequency phase shift, f3 , defined by Eq. (4). The 
transfer function, F, of the normalized, frequency-shifted equivalent 
filter of the ionosphere then becomes, from Eq. (9), 

G(f+f) 2 f 

P(f) = ° = exp(-if o f+ii Z ), J = — (28) 

^ o ; c 

and the associated transfer functions defined by Eq. (15) become 

-if f ? -if f 

A(f) = e — - cos f , A(f) = e — - s ±n f (29) 

When these functions are substituted into the expressions for the 
components of the output spectral density given by Eq. (23), the re- 
sult is 



2 4 2 
9 = f + f 



o o J cp 



W*J(f) = W^C^cos 2 f 2 



Wg(f) = - 4 cos f 2 W^Cf) 



2 

dp W (p) sin 2p 



dp W cp (p)cos(£ 2 +2£ 2 )sin 2 £f 



(f) = 2 J dp W cp (p)W 9 (f-p)sin 2 p(£-f)cos 2 (£ 2 -£f+f 2 ) (30) 
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It is not surprising that the terms in _f appearing in A and A vanish 
in the expressions for the output spectrum, since a term which is 
linear in frequency can be expected to yield only a pure time delay 
to the propagated signal. 

The parameter, f , has the units of frequency, and can be seen 
from Eq. (28) to correspond to that frequency difference from the 
carrier frequency at which the square- law component of the nonlinear 
phase shift produced by the ionosphere equals one radian. Thus, it 
serves as a convenient measure of the degree of nonlinearity in the 
phase characteristic of the ionosphere. Since there is no need to 
consider an accompanying nonlinear amplitude characteristic, f will 
be referred to simply as the !, charac teristic frequency" of the iono- 
sphere at that carrier frequency. 

For computational purposes, it is more convenient to express 
the characteristic frequency in terms of the electron density rather 
than the plasma frequency. As is well known^^ 

2 e 2 N 

^p = T7^ *p = 2nf P < 31 > 

o r 

where 

-19 

e = 1.602 x 10 C is the electron charge 
-31 

m = 9.109 x 10 kg is the electron mass 
12 

€^ = 8 # 854 x 10 F/m is the permittivity of free space 

3 

N = electron density in electrons/m 
e J 

Substituting in Eq. (26) and defining 

N T = J N e ds, electrons/m 2 (32) 
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as the electron density integrated along the propagation path then 
yields 

f 3/2 

f = 1.088 x 10 3 -~ , Hz (33) 

where is the carrier frequency in Hz. The characteristic fre- 
quency is plotted as a function of carrier frequency for a number of 
integrated electron densities in Fig. 1. 
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Fig.l — Characteristic frequency of the ionosphere 



14 



V, COMPUTATIONS FOR FM WITH A UNIFORM BASEBAND 

One transmitted signal of practical interest is that generated 
by using a gaussian baseband signal, with a spectral density uniform 
in (-B,B) and zero elsewhere, to frequency modulate a carrier at an 
rms frequency deviation, D. Such a baseband signal closely approxi- 
mates a multichannel , frequency- di vis ion-multiplexed telephone 
signal. It also approximates a television video signal, though 
less well, partly because typical video spectra are not uniform 
and partly because the amplitude distribution of a video signal is 
less nearly gaussian than that of a multichannel telephone signal. 

The spectral density of the input phase rate will therefore 
be taken as 

V f)=2 ir> |f|sB (34) 

and zero elsewhere. Then, using Eqs. (21) to convert between the 
spectra of phases and phase rates and substituting Eq. (34) into 
Eqs. (30)| yields the following spectral components of the output 
phase rate 

Other related properties are discussed in the next section. 
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— — cos f , |f| £ B 

D Z /f - 
c 

W l (f) 8n 2 ,~r ^ 



D 4 /f^ B 2 o £ 

wf(f) 2rr 2 f 2 "p dp 



2 P d P 2 2 2 II 

cos _f J — cos(jl +2£ )sin _pf , |f | s 



J 



X sin 2 p(£-|f |)cos 2 Cp 2 -p|f j+f 2 ) , |fj ^ 2B (35) 



where it is understood that the terms vanish in those frequency inter- 
vals in which they are not defined and where underlining indicates 

normalization to f . 

c 

The spectral density of the linear- signal component of the output 
phase rate given by the first of Eqs. (35) is plotted as a function 
of the relative baseband frequency, f/B, for a number of normalized 
baseband frequencies, B/f , in Figs. 2 and 3. A feature of note is 
the progressive drop in the high-frequency content as the baseband 
is widened relative to the characteristic frequency. 

The magnitude of the cross-power spectral density given by the 
second of Eqs. (35) is plotted in a similar fashion in Figs. 4 and 5. 
Unlike the linear- signal component which results from a pure auto- 
correlation and therefore has a real, positive, even spectral density, 

These and all subsequent computations were performed numerically 
using the RAND JOSS computer. JOSS is the trademark and service mark 
of The RAND Corporation for its on-line time- shared computer program 
and services using that program. 
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Fig. 2 — Output linear-signal spectral density 
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Fig. 3 — Output linear-signal spectral density, db 
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Fig A — Magnitude of output cross-power spectral density 
(Note that is negative in this example.) 
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Relative baseband frequency, f/B 



Fig. 5 — Magnitude of output cross-power spectral density, db, 
(Note that W? is negative in this example.) 
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the cross-power stems from a cross-correlation and need not be so 
constrained. In this case, the spectrum is real and even but neg- 
ative, at least for the range considered for the parameters. Other 
than that, the principal distinguishing feature of the cross-power 
spectral density is that it generally increases steadily across the 
baseband . 

The intermodulation spectral density given by the third of Eqs. 
(35) is plotted in Figs. 6 and 7. Since the intermodulation is derived 
from an autocorrelation function, it too is a real, positive even func- 
tion. As is characteristic of a second-order spectral density, it 
extends to twice the highest baseband frequency, falling smoothly to 
zero. Like the cross-power, the intermodulation increases across the 
baseband. Since the linear- signal is roughly constant, it follows 
that the cross-correlation between the intermodulation and the linear- 
signal, as indicated by the cross-power, is approximately constant 
across the baseband. 

Narrow-band approximations to these spectra are readily obtained 
by expanding the circular functions in Eqs. (35) and retaining only 
the leading terms. Then, 




|f j £ B 



c 




c 




(36) 
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These spectra are plotted as dashed lines in Figs. 4 to 7 for the lower 



al., mentioned in the Introduction, agree with the third of Eqs. (36) 
above. As seen in Figs. 4 to 7, the approximation is quite good for 
values of B/ f c up to about 0.5. 

The ratio of the linear- signal to the intermodulation spectral 
densities, i.e., the ratio of the first and third of Eqs. (35), is 
plotted in Fig. 8. As explained in Ref. 6, this ratio approximates the 
signal- to-cross-talk ratio in a given channel of a multichannel, frequency 
division-multiplexed modulating signal as a function of its position with 
in the baseband. The approximation improves as the number of channels 
is increased. The poorest channel (i.e., the one with the lowest signal- 
to-cross- talk ratio) lies at the upper end of the baseband, a fact which 
compounds its difficulties since it is also the channel with the lowest 
signal-to-noise ratio due to thermal noise in the r-f channel. 

Also shown in Fig. 8, as a dashed line, is the approximation avail- 
able by forming the ratio of the first and third of Eqs. (36), i.e., 



To this approximation, the signal- to-cross- talk ratio, SCR, in the 
poorest channel is given by 



values of B/f . The results by Rice, Shaft, Prosin, and Denton, et 




T 2 9 i i ' 

W'(f) f Q f (2B-|f|) 




(37) 



SCR 



min 



wj(f) 




1 



(38) 



which is obtained by evaluating Eq. (37) at |f| = B. 
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Relative baseband frequency, f/B 

Fi g# 8 — Ratio of output linear-signal-to-intermodulation spectral densities, 
db. This ratio approximates the signal-to-cross-talk ratio in a narrow 
channel as a function of its location within the baseband 
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The signal-to-distortion ratio, SDR, is the ratio of the total 
linear-signal to total intermodulation obtained by integrating the 
first and third of Eqs. (35) across the baseband and forming their 
ratio, which is plotted in Fig, 9. Again, an approximation can be 
obtained from Eqs. (36) yielding 

SDR ~ U (39) 
5 D B 

which is plotted as a dashed line in Fig. 9. Both the signal-to- 
cross- talk and the signal-to-distortion ratios are seen to be very 
well approximated by Eqs. (37) and (39) for the range of the baseband- 
to-characteristic frequency ratio, B/f , considered. 
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f c - Characteristic frequency of 
ionosphere 
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Fig. 9 — Output signal-to-distortion ratio in total baseband, db 
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VI, EXAMPLES AND DISCUSSION 

The foregoing results can be applied to the problem of linking a 
terrestrial station to an earth satellite at synchronous altitude. The 
propagation path traverses virtually the entire ionosphere in a direc- 
tion which depends on the relative location of the station and the sat- 
ellite and can range from the zenith to the near-horizon. The integrated 

electron density, N^, in Eq. (32) for a path to the zenith at midlati- 

17 2 

tudes through an average ionosphere amounts to about 10 electrons/m , 

but can vary by a factor of 10 either way depending on the time of day, 

the season, the position in the solar cycle, and the geographical lo- 
(8 9) 

cation. ' Similarly, a path near the horizon can yield a value for 

18 19 2 

N^, of from 10 to 10 electrons/m under adverse conditions. 

As an example of a particular application, consider a heavily 
trafficked point-to-point communication link via a synchronous earth 
satellite. Assume the down-link consists of a wide-band, FM transmis- 
sion at a carrier frequency, f , of 4 GHz having an rms frequency de- 

o 

viation, D, of 20 MHz and using a CCIF standardized, 960-channel base- 
band which extends to a maximum baseband frequency, B, of 4.028 MHz. 

A worst-case ionospheric path having an integrated electron density of 
19 2 

10 electrons/m then yields a characteristic frequency, f^, of 87 MHz 

from Fig. 1 or Eq. (33). The normalized maximum baseband frequency, JB, 

becomes 0.0463, and the normalized rms frequency deviation, D, becomes 

0.230, so from Eq. (38), the signal-to-cross-talk ratio, SCR . , in the 

mm 

worst channel becomes 8830 (39.5 db) . 

This example indicates that ionospheric dispersion will not cause 
significant intermodulation distortion in a typical wide-band, FM, 
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communication- satellite link except under only the most severe, and, 

therefore, generally infrequent, conditions. The direct dependence 

of SCR . on the various communication parameters is shown by com- 
min r J 

bining Eqs. (33) and (38), 

12 f6 

SCR m . n = 1.401 x 10 _2_ (40) 

D B N T 

The sixth-power dependence on carrier frequency clearly dominates. 

For example, f need only be halved (-18 db) to virtually overcome 

the effect of a tenfold decrease in N T (-20 db) . For a given baseband 

signal and r-f bandwidth, operation over a given path should be at as 

high a carrier frequency as possible. 

Another example is the use of a satellite to relay TV programs (in 

color or monochrome) either directly to a viewer or to an intermediate 

station for conventional rebroadcast. Again, assume wide-band FM, but 

now with a carrier frequency, f , of 1 GHz and an rms frequency deviation, 

D, of 8 MHz (corresponding to an r-f channel about 50 MHz in width). The 

highest baseband frequency, B, for the 525- line video used in the United 

18 

States is 4.5 MHz . Let the integrated electron density, N , be 4 x 10 
2 

electrons/m on the basis that transmissions in such a broadcast or re- 
lay mode will be confined to well-defined geographical locations and 
therefore be well above the horizon. In this case, the characteristic 
frequency, f , 17.2 MHz and the normalized parameters, B and D, are, 
respectively, 0.262 and 0.465. Then, from Eq. (39), the output signal- 
to-distortion ratio, SDR, becomes 162 (22.1 db) . Alternatively, from 
Fig. 9, SDR = 15. 5 4- 2 x 3. 32 = 22. 1 db, since D = 0.465 (-3. 32 db) . 
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Unfortunately, it is difficult to equate the subjective effect of 
this signal-to-distortion ratio in the output of an ionospherically 
dispersed KM transmission with that of a similar signal-to-noise ratio 
in the output of a conventional vestigial- sideband AM transmission cor- 
rupted by thermal noise. For one thing, the spectral density of thermal 
noise in the video is uniform for the AM transmission, whereas that of 
the intermodulation in the FM transmission goes roughly as the square 
of the baseband frequency. Thus, one might expect a greater effect on 
the high-frequency video content than on the low frequency video content. 
Also, thermal noise is generally uncorrelated with the signal whereas 
the intermodulation can display a strong correlation as evidenced by the 
magnitude of the cross-power spectral density. Therefore, the intermodu- 
lation can be expected to differ in appearance and subjective effect 
from the customary "snow." 

Other aspects, unique to a video signal, must also be considered. 
One is the possibility that the intermodulation may affect the sound sub- 
carrier, which is located above the video at 4.5 MHz. Although the inter- 
modulation spectrum is at its strongest there, the sound subcarrier is 
relatively strong (a minimum of 50 percent of the peak picture power) 
and uses FM with a peak frequency deviation of only 25 kHz. These factors, 
plus the known superior coverage of the sound portion of conventional TV 
broadcast compared with the picture, suggest that this source of inter- 
ference will probably not be significant. 

The situation is less promising with respect to color TV. The 
color subcarrier is at 3.58 MHz with its attendant sidebands extending 
from roughly 2.5 to 4.5 MHz, again in the vicinity of the peak of the 
intermodulation spectral density. As a result, two effects can be 
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anticipated. First, the ratio of the linear- signal- to-intermodulation 
spectral density at which the color information must be demodulated 
(corresponding to the signal- to-cross*- talk ratio in the previous 
example) will be unfavorable, amounting to about 18 or 19 db for this 
case. As a result, the color quality will be degraded much as it 
would be by thermal noise of comparable magnitude. 

Second, the intermodulation can be expected to interfere with 
the synchronization of the color subcarrier generated in the receiver. 
The color subcarrier is suppressed at the transmitter and only short 
bursts (8 cycles of 3.58 MHz) are sent at the beginning of each hori- 
zontal line (i.e., at a 15,750 Hz rate). Because quadrature modula- 
tion is employed, the reference oscillator must be synchronized 
accurately in phase as well as frequency. Thus, intermodulation can 
degrade color separation by causing synchronizing phase errors. 

It is likely that the subjective effect of ionospheric dispersion 
on a multichannel telephone trunk, as discussed above , can be estimated 
with a fair degree of accuracy using the numerical values computed by 
the theory presented here. However, such is probably not the case with 
respect to TV, particularly with color TV. The computations yield 
numbers, which, if interpreted as corresponding to interference produced 
by thermal noise, would be cause for concern. Whether the interference 
caused by ionospheric dispersion will be more or less objectionable in 
comparison is presently a matter of speculation and will probably remain 
so unti 1 experimental tests are conducted. 



31 



Appendix 



SPECTRAL ANALYSIS 



Let the output phase function given by Eq. (18) be written 



e(t) = p + §+f 2 +f 3 +... 



(A-l) 



as suggested in Section III. Substituting into Eq. (19) for Wq and ex- 
panding the result yields the variety of pairings of the $Xf form de- 
fined by Eq. (22). These will be examined in turn to obtain the various 
components of the output phase spectral density, W ft . 



where the expected value is a constant and the Fourier transform of a 
constant is a delta function at zero frequency, i.e., a dc term. The 
coefficient is given as the leading dc term in the first of Eqs. (23). 

$x§ 

This pairing can also be written directly by noting from Eq. (17) 
that $^ is the response when cp is applied to a filter having an impulse 
response, y . The associated steady-state transfer is seen, from Eqs. 
(13) , to be A so 



By definition, this is the linear- signal component of W Q and is the 
second of Eqs. (23). 



3X3 



This first pairing immediately yields 




(A- 2) 



§X<S = |A(f) | 2 W (f) 



(A- 3) 
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pxf 2 

There are two terms of this type in Eq. (19), but they are equiva- 
lent since (3 is not a function of time. Thus, from Eq. (18), 

oo 

pxf = 2(-0 )(-)•£ 3 f Y ,(M-)E[cp(t-u)-f (t)] 2 du (A-4) 
z o z *) i r 

0 

The expected value in Eq. (A-4) has been computed previously and is 

given by Eq. (C-12) of Ref. 10; viz., 



E[cp(t-u)-§ (t)] 2 = f dp W n (p)[jA(p)| 2 -2 A(p)e i2TW + l] (A-5) 



Performing the integration on u, then yields 



pxf = - 2M dp W (p)A(p)A(-p) ( a - 6 > 



'2 o 



9 



where the first and third terms in Eq. (A-5) have disappeared by virtue 
of the second of Eqs. (12). The functions A and A are identified from 
Eqs. (13). 

As with Eq. (A- 2), the Fourier transform yields a delta function 
at zero frequency. Also, the fact that is an even function can be 
used to write 

CO 

0Xf = - (3 6(f) f dp W (p)[A(p)A(-p)+A(-p)A(-p)] 



2 o 
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(3 6(f) [dp W fn (p)[A(p)A(-p)+A*(p)A*(-P)] 



2P 6(f) dp W rn (p)Re A(p)A(-p) (A- 7) 



o 

— ao 



9 
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where the symmetry properties of A and A given by Eqs. (15) have been 
employed. The coefficient of 6 above is given as the second-order dc | 
term in the first of Eqs. (23). 

The computation of the higher order terms requires evaluating 
the expected value of fourth-order products of the gaussian variates. 
Two of the methods available have been used in previous reports. In 
one of these, the desired expected value is identified as the coef~ 
ficient of the appropriate term in the expansion of the expected value 
of a specially constructed exponential function, details are given in 
Appendix B of Ref. 6. The other method uses the relation 



E(x 1 x 2 ...x 2n ) = £ E(x,x;)E(x.x k ) ... E^x^) 
all pairs 

as described in Appendix C of Ref. 10, which is used here. 



(A-8) 



There are two terms of this type in Eq. (19), i.e., 

^Xf^ = 3 E[$(t)f 3 (t+T) + $(t+T)f 3 (t)] 

= 3 E[$(t-T)f 3 (t) + $(t+T)f (t)] (A-9) 

Since the terms in Eq. (A-9) differ only in the sign of t, they will 
simply be written singly with a =t notation, the sum being understood. 
Then, from Eq. (18) , 

CD 

§xf 3 = 3 E {- \ Mt±T) J Y r (x)[cp(t-x) - $ r (t)] 3 dr 

0 

00 

- ^ * (t*T)$.(t) f V, (x)[cp(t-x) - §(t)] 2 dx} (A- 10) 

l r i ti l r j 
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The term in the upper line of Eq. (A- 10) is identical in form to Eq. 
(C-4) of Ref. 10. Its value is given by Eq. (C-18) of Ref. 10 as 

QC 

2 W (f) J dp W^CP) {Re A(p)A(-p-f)A(f) - | A (f) | 2 | A (p ) \ 2 ] (A-ll) 

- CO 

The other term in Eq. (A- 10) remains to be computed. Using Eq. 
(17) to expand the initial s yields 



CO CO CO 

r» 



1 r r ^ 

2 J du J dv J dw \( u )'Y i (v)Y i ( w ) 
0 0 0 



3 E cp(t±T-u)cp(t-v)[cp(t-w) - * r (t)] 2 (A- 12) 

According to Eq. (A-8) , the expected values in this expression can be 
expanded in the form 

E(abc 2 ) = E(ab)E(c 2 ) + 2E(ac)E(bc) (A- 13) 

The various expected values are 

E(ab) = E cp(t±T-u)cp(t-v) 

= R (±t-u+v) (A- 14) 

9 

where the autocorrelation function in Eq. (20) has been used, 
E(c 2 ) = E[cp(t-w) - $ r (t)] 2 

CO 

dp Wcp (p)[|A(p)| 2 - 2A(p)e i2TIPw + 1] (A- 15) 



which follows from Eq, (A-5) , 
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E(ac) = E cp(t±r-u) [cp(t-w) - $ (t)] 

CO 

= R^iT-ufw) - J dx V r (x)^(±rr-vH-x) 



E(bc) = E tp(t-v) [cp(t-w) - 4 (t)] 

00 

= R (v-w) - J dx v r (x)R (v-x) 
0 

= J dp W ( P )e i2n e< V - W > - J dx Y (x) J dp w (p)e 12n P< v 

°£ i2npv r- 

- \ *> W^(P) e [e" i2n P W - A( P) ] 



The Fourier transforms of Eqs. (A-14) and (A-16) are 



3 E(ab) = J dT e" l2nfT R^C^r-ufv) 



E(ac) = J dT e" l2ITfT [^(iT-u+w) - J dx Y r (x) Rp (±T-u+x) J 



Ti2nf(u-w) „ Ti2nf u . 

= e WpCf) - e A(Tf) Wcp (f) 



TiZnfu 

= e 



Combining Eqs. (A- 15), (A- 17) - (A- 19) in Eq. (A- 13) yields 
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3 E(abc 2 ) = E(c 2 )5 E(ab) + 2E(bc)5 E(ac) 



Ti2nfu, 

= e 



+ 2 e 



w (£) { e *2nfv J dp Wcp(p) [| A(p) |2. 2A(p)e i2"Pw +1 ] 

_ 00 

[e ±i2TTfw -A(Tf} J dp W cp (p)e i2TTpv [e- i2TTpw -A(p)]} (A-20) 



When substituted into expression (A-12), the integration on u follows 
immediately , resulting in 

CO CO 

1 * / t /■ \ f i r , / n f \ fterm in braces\ /A 01N 

- 2 A ( ±f > W cp( f ) J dv J dw YiMYiW { in Eq . (A _ 20 ) } < A " 21 > 
0 0 

Integrating on v further yields 

GO 

- -| A(±f)W (f) J dw YjL (w) 
0 

CO 

W) J dp W (p)[|A( P )| 2 -A( P )e i2TTpW + l] 

-CO * 

+ 2 [e ±i2TTfw -A(*f)] J dp W cp ( P )A(-p)[e- i2TTpw -A( P) ] } 



~ CO 

- \ A(±f) Wcp (f) J dp vyp) J dw Y.(w) 

-co 0 

{A(^f)[|A(p)| 2 -2A(p)e i2TTpw + l] 



+ 2A(-p)[e ±i2TTfw -A(Tf)] [e- i2TTpw -A(p)] } (A-22) 

The terms not involving the exponentials vanish when integrated on w, 
by virtue of the second of Eqs. (12), yielding 
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\ A(±f) Wcp (f) J dp typ) {-2A(Tf)A(p)A(-p) 

— OO 

+ 2A(-p)[A( : Ff + p).A( : Ff)A(p)-A(p)A(Tf)]j. 



dp W tp (p)A(±f)[2A(Tf)A(p)A(-p) 



- A(-p)A(Tf+p)+A(-p)A(Tf)A(p)] 
= W^(f) J dp W^CP) {A(±f)A(-P)[2A(Tf)A(p)-A(Tf+p)] 

-00 

+ |A(^f)| 2 |A(p)| 2 } (A-23) 

where the symmetry properties of A and A given by Eqs. (15) have been 
used in the last step. Dropping the ± notation and writing expression 
(A-23) as a sum yields 



V f ) 



dp W rp)[2A(f)A(-p)A(-f)A(p) - A(f)A(-p)A(-f+p) 

T 
- 

+ |A(f) | 2 |A(p) | 2 + 2A(-f)A(-p)A(f)A(p) 

- A(-f)A(-p)A(f+p) + |A(-f) | 2 |A(p) j 2 ] (A-24) 

Changing the sign of the variable, p, in the second and third terms 
results in 

W rf (f) f dp W r (p)[2A(f)A(-p)A(-f)A(p) - A(f)A(p)A(-f-p) 

— QC 

+ 2A(-f)A(p)A(f)A(-p) - A(-f)A(p)A(f+p) + 2 |A(f) | 2 |A(p) | 2 ] 

CO 

= 2W (f) f dp W fp)[2 Re A(f)A(-p)A(-f)A(p) - Re A(f)A(p)A(-f-p) 
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+ 2|A(f)| 2 |A(p)| 2 ] (A-25) 
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where the symmetry properties of A and A have again been used. 

The total $Xf^ i s the sum of expressions (A-ll) and (A-25) , 

CO 

$Xf 3 = 2W (f) J dp lyp) {lie A(f)[A(p)A(-p-f) - A(p)A(-f-p)] 

— CO 

+ |A(f)| 2 [|A(p)| 2 -|A(p)| 2 ] + 2Re A(f)A(~f) A(p)A(-p)} (A-26) 

This contribution is the leading cross-power term and is listed as the 
third of Eqs. (23) 



f 2 Xf 2 



From Eq. (18), this term is 

CO 00 

f 2 Xf 2 = I I du Y i (u) J dv V V) 
0 0 

3 E{[cp(t-u) - $ r (t)] 2 [cp(t+T-v) - ^ r (t+T)] 2 } (A-27) 

Expanding the expected value according to Eq. (A- 8) yields 

E(a 2 b 2 ) = E(a 2 )E(b 2 ) + 2E 2 (ab) (A- 28) 



The expected values 



E(a ) = E[cp(t>u) - $ (t)] 



r 



E(b ) = E[cp(t+T-v) - * (t+T)] 2 (A-29) 

are equivalent and are both given by Eq. (A- 15). Since the variable 
T disappears on applying the expectation operator, taking the Fourier 
transform results in a zero-frequency delta function. Thus, the con- 
tribution of this term to Eq. (A-27) is 

CO co CO c° 

■J 6(f) j du Y .(u) j dv y ± (v) J dp W (P) J da w (a) 

0 0 

[|A(p)| 2 -2A(p)e i2npu +l] [|A(a)| 2 _2A(a)e i2npv + l] (A- 30) 
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Again by virtue of the second of Eqs. (12), only the terms involving 
the appropriate exponentials remain when the integrations on u and v 
are performed. The result is 

CO 00 

6(f)[ J dp Wcp (P)A(p)A(-p)] [ J da W (a)A(a)A(-a)] 

-GO — OO 

oo 2 

= 6(f)[ J dp W (p)Re A(p)A(-p)l (A-31) 

-00 

where the symmetry properties of A and A are again used. This term 
is not listed in Eqs. (23), since it is clearly a third-order dc term 
(another contribution of the same order can be expected from the 
^ xf term) . 

The expected value in the second term of Eq. (A-28) is similar 
in form to Eq. (C-26) of Ref. 10, 

00 

E(ab) = J dp W cp (P)e i2npT [A(p)-e- i2rTpV ] [a ( -p)-e i2rTpu ] (A-32) 

— OO 

Substituting first into Eq. (A-28) and then into Eq. (A-27) yields 

00 00 

2 J du I dv ^i( u ^i( v > 

CO 
—CO 

r Arr t t ,„\ i^rA/.N -i2nav*l [\ , ^ i2nau! \ 
J da w (a)e A(a)-e A(-a)-e J 

—00 ^ 

-iJdpJdc„<p) V a) 3 e 12 ^«> T 

-OO 

{Jdu Yi (u)[A(-p).e i2npu ][A(-a)-e i2nCTu ] 

J dv Yi (v)[A(P)-e- i2npV ] [A(a)-e- i2nCTV ] } (A- 33) 



0 0 



0 

00 
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The integral on u in the second line is evaluated by expanding the 
integrand and integrating term by term 

J du Y .(u) LA(-p)A(-a)-A(-^ 
0 

= - A(-a)A(-p)-A(-p)A(-a)+A(-p-a) (A- 34) 

Analogously, the integral on v in the third line is 

-A(a)A(p) - A(p)A(a) + A(p+^) (A-35) 

The product of these expressions, which are complex conjugates, is 
simply the square of their magnitude. Substituting in (A-33), 
taking the Fourier transform and integrating on o yields 

^ J do J da W.^(P)fcL(c)6 (f-p-a) 

_CO _CO ^ ^ 

|A(p)A(a) + A(p)A(o) - A(p^r)| 2 
" \ J dp W cp (P)W cp (f-P)|A(p)A(f-p)+A(p)A(f-p)-A(f)| 2 (A-36) 

— co 

which is the leading intermodulation component and is listed as the 
last of Eqs. (23) . 
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